Polyakov Loop Models, Z(N) Symmetry, and Sine-Law Scaling 
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We construct an effective action for Polyakov loops using the eigenvalues of the 
Polyakov loops as the fundamental variables. We assume Z{N) symmetry in the 
confined phase, a finite difference in energy densities between the confined and de- 
confined phases as T ^ 0, and a smooth connection to perturbation theory for 
large T. The low-temperature phase consists of — 1 independent fields fluctu- 
ating around an explicitly Z{N) symmetric background. In the low-temperature 
phase, the effective action yields non-zero string tensions for all representations with 
non-trivial A^-ality. Mixing occurs naturally between representations of the same 
A^-ality. Sine-law scaling emerges as a special case, associated with nearest-neighbor 
interactions between Polyakov loop eigenvalues. 

Effective actions for the Polyakov loop are directly relevant to the phase diagram and 
equation of state for QCD and related theories [l|, There is also good reason 

io believe that Polyakov loop effects play an important role in chiral symmetry restoration 

^ l^- Furthermore, there may be a close relatioriship between the correct effective 
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action and the underlying mechanisms of confinement jlll . 

In this letter, we construct a general effective action for the Polyakov loop, making as 
few assumptions as possible. Our goal is a comprehensive model which describes both 



the confined and deconfined phases of SU{N) pure gauge theories at finite temperature, 
as well as the properties of the deconfining phase transition. We will describe below the 
general form of such an effective action, using the Polyakov loop eigenvalues as natural 



variables. Effective actions of this type have previous 



temperature 



strong coupling 
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KJP in two dimensions 



been constructed for QCD at high 
2Q|, and for lattice gauge theories at 



The Polyakov loop P{x) is the natural order parameter for the deconfinement transition. 
It is defined as a path-ordered exponential P{x) = Pexp i (ItAq {t,x) where xis a d- 
dimensional spatial vector and (3 is the inverse temperature. The confined phase has a Z[N) 
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symmety that implies {Trp (a;)y = for all n not divisible by A^: k\N ^ For these 
powers of the Wilson line, the asymptotic behavior 

TrpP^ (x) TrpP^^ {y)) oc exp [-(3ak \x - y\] (1) 

is observed in lattice simulatiosn as |x — ?/| — » oo. The string tension ak is in general 
temperature-dependent. Similar behavior is also seen for the two point functions associated 
with irreducible representations of the gauge group {Tr^P (x) Tr^P^ (y)). Every irreducible 
representation R has an associated A^-ality kji such that Tr^iP — > z^^Tr^P under a global 
Z{N) transformation, with G {0,..,A^ — 1} and z G Z{N). Any representation with 
kn ^ gives an operator Tr^P which is an order parameter for the spontaneous breaking 
of the global Z{N) symmetry. 

It is widely held that the asymptotic string tension depends only on the A^-ality of the 



representation. In fact, lattice data 



do not yet show this behavior 
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rom simulations of four dimension SU (3) gauge theory 



26|. Instead, the string tension an associated with a 



representation R scales approximately as 

<yR = yr^^F (2) 

where Cr is the quadratic Casimir invariant for the representation R. This behavior can be 
rigourously demonstrated in two-dimensional gauge theories. For simplicity, we will refer to 
this as Casimir scaling. We use the term Z{N) scaling to descibe the asymptotic behavior 

k{N-k) 

an = i^-. (3) 

where k is the A^-ality of the representation R. As we will show below, Z{N) scaling is 
obtained from Casimir scaling at large distances if there is mixing between representations 
of the same A^-ality. If Z{N) scaling holds, there are [A^/2] distinct asymptotic string 
tensions, where [A^/2] is the largest integer less than or equal to A^/2. For the first [A^/2] 
antisymmetric representations made by stacking boxes in Young tableux, Casimir and Z{N) 
scaling are identical. Another possible scaling law for the string tensions is sine-law scaling 

sm(^ 



which has been shown to occur in softly broken N = 2 super Yang-Mills theories|27] and in 



MQCD 
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In a gauge in which Aq is time independent and diagonal, we may write P in the funda- 
mental representation as 

Pjk = exp {iOj) 6jk (5) 

where we shall refer to the numbers 6j as the eigenvalues. They are not independent 
because det (P) = 1 implies 

Y,0j = O mod 2tt. (6) 

j 

The information in the different representations is redundant. All the information is con- 
tained in the — 1 independent eigenvalues of P. 

Without spatial gauge fields, the only gauge-invariant operators we can construct are class 
functions of P, depending solely on its eigenvalues. Thus the effective action should take 
P to be in the Cartan, or maximally commuting, subgroup U{1)^~^. An effective action 
constructed from SU{N) matrices would introduce spurious Goldstone bosons associated 
with the off-diagonal components of P in the deconfined phase. On the other hand, the 
simplest Landau-Ginsburg models of the SU (2) and SU (3) deconfining transitions have used 
TrpP as the basic field. These two cases are special because TrpP specifies the Polyakov 
loops in all other representations. In SU{4), there are sets of eigenvalues for which 

TrpP = and TrpP^ = (7) 

consistent with a Z{A) symmetry, and another set for which 

TrpP = and TrpP^ ^ (8) 

consistent with a Z{2) symmetry. Thus knowledge of TrpP alone is insufficient for > 3^1. 
From the characteristic polynomial, one may show that the eigenvalues of a special unitary 
matrix are determined by the set ^TrpP^^ with k = 1..N — 1, and of course vice versa. 

In both the confined and deconfined phases, we would like our effective theory to proceed 
from a classical field configuration which has the symmetries of the phase. If we denote that 
field cofiguration as Pq in the confined phase, Z{N) symmetry requires that TrpP^ = for 
all k not divisible by A^. Enforcing this requirement for = 1 to A^ — 1 leads to a unique set 
of eigenvalues via the characteristic equation z'^ + (—1)^ = 0, but it is instructive to derive 
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the set another way. For temperatures T below the deconfinement transition T^, center 
symmetry is unbroken. Unbroken center symmetry imphes that zPq is equivalent to Pq after 
an SU (N) transformation: 



zPo = gPog+. (9) 

This condition in turn implies TtjiPq = for all representations R with non-zero A^-ality, 
which means that all representations with non-zero A^-ality are confined. The most general 
form for Pq may be given as hdh^, where h G SU{N), and d is the diagonal element of 
SU{N) of the form 

d = w diag z,z^,..,z^ = 1 (10) 

where z is henceforth exp (27Ti/N), the generator of Z{N). The phase w ensures that d has 
determinant 1, and is given hj w = exp [— (A^ + 1) ni/N]. Strictly speaking, w is required 
only for N even, but it is convenient to use it consistently. We will henceforth identify Pq 



with d. Another useful representation is {Po)jk = ^jfcexp 



where 



TT 



^? == ^ (2j - AT - 1) . (11) 

Thus the Z(A^)-symmetric arrangement of eigenvalues is uniform spacing around the unit 
circle. This is consistent with the known large- A^ behavior of soluble models in the confined 



phase 



29|. 



Assume that Pq is the global minimum of the potential V associated with the effective 
action for temperatures less than the deconfining temperature T^. Because the gauge fields 
transform as the adjoint representation, is a class function depending only on represen- 
tations of zero A^-ality. It is thus a function only of the differences in eigenvalues 6j — 6^, 
with complete permutation symmetry as well. In the low-temperature, confining phase, we 
consider small fluctuations about Pq, defining 6j = 6^ + 66 j. Although this approximation 
may not be a priori valid, the assumption that fluctuations are small can be justified in the 
large- A^ limit. For small fluctuations 

N N 

TrpP^ = w^z^^'e'^^^- ~ ^ w^z^'^ik 56.^ (12) 

n=l n=l 

provided k is not divisible by A^. Note that TrpP^ takes the form of a discrete Fourier 
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transform in eigenvalue space. We define the Fourier transform of the fields 0„ as 

N 
n=l 



(13) 



so TrpP^ = ikw^(j)k- We will show below that the fields are the normal modes of the 
Z(A^)-symmetric phase in a quadratic approximation. The mode (pN = <Po is identically zero 
for SU{N), and from the reality of 6, we have (pN-n = (0n)*- In the case where k is divisible 
by A^, TrpP^ has a leading constant behavior of Nw'', and the term linear in (p vanishes. 
The adjoint Polyakov loop operator is approximately Tr^P = TrpP TrpP^ — 1 ~ 010^ — 1. 
The operator (picpX has a non-zero vacuum expectation value, and should couple to scalar 
glueball states. 

Operators with the same A^-ality generally give inequivalent expressions when written 
in terms of the the variables. For example, TrpP^ oc 0^ and [TrpP)^ oc (0i)'^. Group 
characters Xr{P) are represented as sums of terms with the same A^-ality, but with different 
mode content. For example, in 5'[/(4), the 10 and 6 representations are given by 



XS,A — - 



{TrpPf ±TrpP^' 



±202 + i^l 



(14) 



As we discuss below in detail, these different combinations of fields will in general produce 
several different excitations, and only the lightest states will dominate at large distances. 

The form of the effective action is fixed at high temperature by perturbation theory. The 
form of the effective action at high temperatures can be written as 
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S,ff = (3 j (fx [T'^Trp {Vef + ¥^^{9) 



(15) 



The kinetic term is obtained from the underlying gauge action via ^TrpF 
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\2Trp (VAc 



T'^Trp (V6')^. The potential Vil{6) is obtained from one-loop pertur- 



bation theory. For our purposes, it is conveniently expressed as 
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n=l 



iTrP'^r - 1 



g 2 



n=l 



7r2 n4 



AT-l + ^cos {n (6,-6^)) 



This series can be summed to a closed form in terms of the 4th Bernoulli polynomial. The 
complete one- loop expression has been obtained recently 3^0, 34 1; the complete kinetic 



term has a 6'-dependent factor in front of the derivatives. 



6 



There are N equivalent solutions of the form 



er = -jf (16) 

related by Z{N) symmetry breaking. All of these solutions break Z{N) symmetry, with 
TrpP — N exp {2'Kip/N). For these values of 6, we recover the standard black-body result 
for the free energy. 

A sufficiently general form of the action at all temperatures has the form 



s^ff ^(5 j Sx [kt'^ttf {vef + v{e)\ (i?) 



where « is a temperature dependent correction to the kinetic term, and y is a function 
only of the adjoint eigenvalues 9j — 9k- More complicated derivative terms can be added as 
necessary. 

We assume that there is a finite free energy density difference associated with different 
values of P as T ^ 0. Because the eigenvalues are dimensionless, this requires terms in the 
potential with coefficients proportional to {mas s)'^ as T — > 0. We can expand the potential 
to quadratic order around Pq 



v{e)c.v{e') + j:l 



dOjdek 



d9jd9k (18) 



where the coefficient in the expansion depends only on \j — k\. The quadratic piece is thus 
diagonalized by Fourier transform, and we can write 

N-l 

V {9) ~ V{9') + ^ M^cPn'pN-n. (19) 

n=l 

Similarly, the kinetic term becomes 

j2 N-l 

kT^Ttf {V9f = ^ E (V*^") (V0iv-n) . (20) 

n=l 

Once an ordering of eigenvalues is chosen, Z{N) symmetry is expressed as a discrete trans- 
lation symmetry in eigenvalue space. If we write the higher-order parts of Se/f in terms 
of the Fourier modes each interaction will respect global conservation of A^-ality. For 
example, in 5'C/(4), an interaction of the form 0^02 is allowed, but not 0i02- 

The confining behavior of Polyakov loop two-point functions at low temperatures, for 
all representations of non-zero A'"-ality, is natural in the effective model. If the interactions 
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are neglected, we can calculate the behavior of Polyakov loop two-point functions at low 
temperatures from the quadratic part of S^ff. We have for large distances 



where (T„ (T) = y NM^(T) / k(T) is identified as the string tension for the n'th mode at tem- 
perature T. Interactions may cause the physical string tensions to be significantly different 
from the tree- level result, but the general field-theoretic framework remains in any case. Of 
course, 4>N-n = {(pn)* implies an (T) = crAr_„(T). The number of different string tensions is 
[A^/2], the greatest integer less than or equal to N/2. The zero-temperature string tension 
is given at tree level by 



This formalism gives a natural mechanism for the transition from Casimir scaling to scal- 
ing based on A^-ality at large distances. As we have noted above, there are many composite 
operators with the same A^-ality. For example, in SU{8), the operators (pf, 0f02, 01035 
and 04 all have A^-ality 4. Naively, the string tensions associated are 4cti, 2ai + a2, 2a2, 
c"i + c"3; and (74, respectively. However, there is no symmetry principle prohibiting mixing of 
operators of the same A^-ality. If there is mixing, then only the lightest string tension will 
be observed at large distances. Interactions between modes will lead to such mixing, and 
the operators TrpP'^ will exhibit a more complicated behavior. 

Mixing between Polyakov loop operators of the same A^-ality is easily understood using 
character expansion techniques applied to lattice gauge theories. For d > 2, there are strong- 
coupling diagrams in which the sheet between two Polyakov loops in a representation R split 
into a bubble with sheets of representation Ri and i?2- Such strong-coupling graphs are non- 
zero when i? C -Ri ® -R2 ; a necessary but not sufficient condition is that R and Ri ® R2 have 
the same A^-ality. The same graph couples two representations R and R' of the same A^-ality 
via the process R Ri + R2 R' ■ The case ci = 2 is special: in the continuum, there are 
no gauge vector boson degrees of freedom, and the theory can be solved exactly, yielding 
Casimir scaling. The lattice theory is also exactly solvable; it reduces to a = 1 spin chain 
model, and there are no bubble diagrams as in higher dimensions. Although the string 
tensions associated with a given lattice action do not in general give Casimir scaling, the 
fixed-point lattice action in two dimensions does yield results identical to the continuumj^^. 



(Tr^P" {x) TrFP+^ {y)) cx (0„ (x) 0^ {y)) cx exp 



(21) 




al (0) = ArM„^(0)/«:(0) 



(22) 
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It is easy to see that the [A^/2] string tensions ai, (72, .., cr[Ar/2] are all set independently 
within the class of effective models. A minimal model for the confined phase exhibiting 
this behavior is 

[N/2] ^4 

V =Y, -^TrFP^TrpP^^ (23) 
k=i ^ 

where the are arbitrary. The /c'th term in the sum forces TrpP^ = 0, and gives rise to 
a mass for the mode (pk- 

Sine-law scaling arises naturally from a nearest-neighbor interaction in the space of 
Polyakov loop eigenvalues. Consider the class of potentials with pairwise interactions be- 
tween the eigenvalues 

V2 = j:Hd,-0k). (24) 



as obtained, for example, by two-loop perturbation theory 
shows that at tree level 



An elementary calculation 



(Tr, 



where v^"^^ is the second derivative of v. This master formula relates the string tensions to the 
underlying potential. It is essentially the dispersion relation for a linear chain with arbitrary 
translation-invariant quadratic couplings: nearest-neighbor, next-nearest neighbor, et cetera. 
If the sum is dominated by the j = 1 and j = N — 1 terms, representing a nearest-neighbor 
interaction in the space of eigenvalues, then we recover sine-law scaling 



m Sin m . (26) 



There is a large class of potential which will give this behavior. Z{N) scaling can be obtained 
by a very small admixture of other components of v^'^\ as we show below. 



3een measured in d = 3 and 4 
We examine the 



The string tension associated with different A^-alities has 

dimensions for = 4 and 6, and in ci = 4 for A^ = 8 0, y, H Q, Q 

simulation data by inverting equation ()25|) above to give a measure of the relative strength 
of the couplings v^'^^ {2Tcj /N). We normalize the result of this inversion such that the sum of 
the independent couplings adds to one, and the results are shown in Table HJ Sine-law scaling 
corresponds to a value of Ifor j = 1, and for the other [A^/2] — 1 independent couplings. 
For Z{N) scaling, the large- A^ limit gives f*^^-* {2tx j/N) oc yielding the result shown in 
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TABLE I: Relative strength of couplings w^^^ {2TTj/N) 







7 = 1 


7 = 2 


7 = 3 

J 


7 = 4 


kJ y-t 1 \kj 1 J (JO 


KJ 




043l'9'l 






QTr(A)\A^ 1 


A 


n Qfisnfi"! 


0S9('7') 






'^TKdM'^Q 4nl 






onsi'^') 






'^TI(d) 7,(N\ 




n Q41 9 


0^88 






kJ Ly \\J J \kj 1 J OO 


(J 


n Q'^Di'^'i 


06^l'8') 


n no4('5'i 






A 


n QfiDi'ifi'i 


04^('1 R\ 










n QQfiMn'i 


0n0'^('91 R\ 


n nodOP,') 




SU(6) Z(N) 


any 


n 9266 


0.0618 


0.0116 




5C/(8)f4]J 


4 


1.028(22) 


-0.067(24) 


0.047(32) 


-0.009(22) 


5C/(8) Z(iV) 


any 


0.9249 


0.0583 


0.0130 


0.0037 


Z{N) N ^oo 


any 


0.9239 


0.0577 


0.0114 


0.0036 


sine Law 


any 


1 












the table. Note that the difference between sine-law scaling and Z{N) scaling remains small 
but finite, even as goes to infinity. The three-dimensional simulation results clearly favor 
Z{N) scaling. The two sets of four-dimensional simulation results for SU{A) and SU{6) 
agree within errors, but one set lies systematically closer to the sine-law predictions, while 
the other set of results does not really favor either theoretical prediction. The SU (8) results 
nominally favor the sine-law prediction. However, smaller error bars will be necessary to 
differentiate between sine-law and Z{N) scaling in four dimensions. 

The ultimate origin of confinement influences the form of the potential V, and the be- 
havior of the different string tensions are derived in turn from V. In previous work on 
phenomenological models of the gluon equation of statej^ Is], we considered models of the 
form 



f = -p= V{9) - 2 



(2^ 



Tta In 



1 - Pe-'^'^^ 



(27) 



where is a phenomenologically chosen potential whose role is to favor confinement at 
low temperature. We studied two physically motivated potentials which reproduce SU{3) 
thermodynamics well. Both give a second-order deconfining transition for SU{2), and a first- 
order transition for SU {N) with > 3 in accord with simulation results. One potential is 
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a quadratic function of the eigenvalues 

N a~l 

Va = vaJ2T. i^o. - Op) {9^ -Op- 27r) (28) 

o=2 /3=1 

which appears as an 0{m?T'^) term in the high temperature expansions at one-loop. This 
potential leads to = ai for every A^-ality. The other potential is the logarithm of Haar 
measure 

N a-1 

= E E In [1 - cos {9^ - 9p)] . (29) 

a=2 13=1 

which is motivated by the appearance of Haar measure in the functional integral. This term 
is cancelled out in perturbation theory for fiat space 17 1, but not in other geometries [2^. The 



potential Vb was first studied by Dyson in his fundamental work; on random matrices 
and leads to 

o I kjN-k) 
= V iV-1 

which one might call "square root of Z(N) scaling". Although these two models are not 
consistent with the lattice simulation data for ak for > 3, they remain viable phenomeno- 
logical forms for N = 2 and 3. It is interesting to note that there is a well-studied potential 
which gives Z{N) scaling 43|. It is the integrable Calogero-Sutherland-Moser potential 

(^) = E -ttW (31) 



which has been associated with two-dimensional gauge theoryjl 

The effective action has domain wall solutions in the deconfined phase, generalizing the 
behavior seen in the perturbative, high-temperature form of the effective action |30l. l31|. In 
this limit, the Z{N) symmetry brealcs spontaneously, with N equivalent vacua characterized 
by P = 2;"/. There is a surface tension associated with one-dimensional kink solutions of 
the effective equations of motion which interpolate between different phases. There are 
different surface tensions pk, each associated with the kink solution connecting the n = 
phase with the n = have given an argument 

valid at high temperature indicating that the surface tensions pk obey Z{N) scaling. This 
argument can be extended with minor modifications to arbitrary potentials of the form V2, 
giving semiclassically 

= N-l (^^^ 



11 



in the entire deconfined phase. 

A similar argument has been given for the string tension in the confined phase of 2 + 1- 
dimensional Polyakov models, which consist of gauge fields coupled to scalars in the adjoint 
representation j4^. In this case, the effective potential has the form V2 when written in terms 
of dual variables. It is natural to speculate that there is a class of self-dual effective models 
in two spatial dimensions with identical string-tension scaling laws in both phases. 

We have constructed the most general class of effective actions which can be built from 
the eigenvalues of the Polyakov loop. Z{N) symmetry requires a symmetric distribution 
of eigenvalues in the confined phase, leading to a special role for Z{N) Fourier modes. 
Mode mixing provides a natural crossover from Casimir scaling to scaling based on A^- 
ality. There is a natural association between sine-law scaling and strong nearest-neighbor 
interactions between eigenvalues, and sine-law scaling and Z{N) scaling are in some sense 
quite close. Nevertheless, string tensions in different A^-ality sectors are determined by the 
parameters of the effective action, and are a priori arbitrary. It is clear from studies of 
different lattice gauge actions in two dimensions that the ratios of lattice string tensions 
need not be universal. Only near the continuum limit are two-dimensional string tension 
ratios universal. It remains mysterious why the wide class of models so far studied have not 
provided us with a wider variety of string tension scaling laws. We do not know if additional 
fields of zero N-ality can change string tension scaling in the continuum, and there has not 
been much exploration of the possible effect of fields that preserve a non-trivial subgroup of 
Z{N). Studies of such models might give insight into the connection between confinement 
mechanisms and string-tension scaling laws. 
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